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Abstract 

In a first part the scope of classical thermodynamics and statistical mechanics is 
discussed in the broader context of formal dynamical systems, including computer 
programmes. In this context classical thermodynamics appears as a particular the- 
ory suited to a subset of all dynamical systems. A statistical mechanics similar to 
the one derived with the microcanonical ensemble emerges from dynamical systems 
provided it contains, 1) a finite non-integrable part of its phase space which is, 2) er- 
godic at a satisfactory degree after a finite time. The integrable part of phase space 
provides the constraints that shape the particular system macroscopical properties, 
and the chaotic part provides well behaved statistical properties over a relevant 
finite time. More generic semi-ergodic systems lead to intermittent behaviour, thus 
may be unsuited for a statistical description of steady states. 

Following these lines of thought, in a second part non-extensive hierarchical sys- 
tems with statistical scale-invariance and power law interactions are explored. Only 
the virial constraint, consistent with their micro dynamics, is included. No assump- 
tions of classical thermodynamics are used, in particular extensivity and local ho- 
mogeneity. In the limit of a large hierarchical range new constraints emerge in some 
conditions that depend on the interaction law range. In particular for the gravita- 
tional velocity-site scaling relation is derived which is consistant with the 
ones empirically observed in the fractal interstellar medium. 
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1 The scope of classical statistical physics 

1.1 The non-universality of classical statistical mechanics 

The thermodynamical principles were sometimes believed to be applicable to 
all "sufficiently large" physical systems, including the whole Universe. The 
principles of thermodynamics have been so well verified in terrestrial condi- 
tions that the extrapolation of universal validity was adopted by most scien- 
tists. 

In fact this optimism had to be seen as exaggerated after, among others, the 
pioneer works in astronomy of Henon (1961), Antonov (1962), and Lynden- 
Bell & Wood (1968) about the thermodynamics of gravitating systems. These 
works showed that well established rules, such that specific heat is always 
positive, are not verified by all large systems. The awareness that the scope 
of classical thermodynamics is restricted to a subset of physical systems has 
spread slowly but steadily among physicists and chemists (Lynden-Bell 1999). 
In particular a large part of classical thermodynamics is restricted to extensive 
systems. Thus the old assumption, made among others by Kelvin, that ther- 
modynamics in the present form applies to the whole Universe remains to be 
verified. In practise astronomers know well that classical statistical mechanics 
does not work for gravitational systems, since they need to perform intensive 
computer simulations of gravitating particles to predict their average collec- 
tive behaviour and evolution. In contrast, chemists and physicists can rely on 
thermodynamics to predict a large range of possible collective states adopted 
by molecules. 

Therefore it seems useful today to better characterise the effective scope of 
classical statistical physics, and eventually to try to extend it to more systems, 
including non-extensive systems and non physical systems. Such extensions 
toward non-physical systems have been attempted many times, for instance 
in information theory (Shannon 1948), or for image processing. But the as- 
sumptions made are not better founded than the ones of statistical physics, 
therefore much remains to be clarified. 

On the other hand statistical tools are used in many different areas of hu- 
man activities to describe complex systems, natural inorganic systems as well 
as organic or sociological systems. Besides general mathematical tools there 
does not seem to exist uniform statistical principles similar to the ones of 
thermodynamics that can be applied for all the possible systems. 

We will suggest in the following that the main characteristics that systems 
need for applying statistical tools similar to thermodynamics is a chaotic mi- 
croscopic dynamics leading to a fast ergodic behaviour, which may however 
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still be subject to global constraints associated with conservation laws. Suf- 
ficiently chaotic dynamical systems that nevertheless are constraint by sym- 
metry invariants explore a subset of their phase space in a sufficiently good 
ergodic fashion. This allows both to recognise general symmetries coming from 
the phase space restriction, and to use statistics appropriately over the almost 
ergodic domain. 



1.2 Thermodynamics is not necessarily fundamental 

The principles of thermodynamics were elaborated during the 19 th century 
and found empirically to be useful for particular systems, such as thermal 
machines, or diluted gases. Later during the 20 th century the thermodynamical 
principles could be shown to be consistent with a microscopic description 
based on classical or quantum mechanics, provided some assumptions (ergodic 
principle, extensivity, etc.) were introduced. In textbooks these assumptions 
are often not all explicitly stated, in particular that thermodynamical variables 
must be either intensive or extensive. This later rule has been sometimes 
viewed as an additional principle of thermodynamics. 

But some systems, such as turbulent flows, semi-chaotic mechanical systems, 
or systems in a phase transition do not follow globally thermodynamics; yet 
they display reproducible behaviours in a statistical sense. Such systems do 
not belong to systems describable with the tools of statistical mechanics due 
to assumptions being violated. For example, some systems are out of mechan- 
ical or chemical equilibrium, some others are non-extensive, some others are 
insufficiently chaotic. 

So the scope of thermodynamics must be viewed as restricted to particular 
physical systems. The key point to observe is that in certain conditions it is 
possible to simulate thermodynamics on computers, that is, on formal systems. 
This shows that thermodynamics may emerge from purely abstract dynamical 
systems, so is not more fundamental than the underlying dynamical system. 
Molecular dynamics using classical dynamics is today an important formal 
tool to understand natural processes involving millions of atoms. The pioneer 
work by Fermi, Pasta and Ulam (1955) failing to reproduce on computer the 
expected approach to thermal equilibrium by a chain of coupled non-linear 
oscillators did precisely show that the chosen mathematical model may also 
be in conflict with classical thermodynamics. So one can sometimes simulate 
thermodynamics with purely numerical models based on ordinary differential 
equations, which are based for example on classical mechanics; or, it may 
happen that particular classical mechanical systems violate the usual thermo- 
dynamical laws. Thus it seems clear today that classical statistical physics can 
hardly be supposed to be more fundamental than the underlying microdynam- 
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ics of the considered systems. For example the statistical behaviour of a stellar 
cluster has nothing to do with the physical nature of its stars, since the same 
behaviour can be simulated with the classical dynamics of point particles on a 
computer. Instead of a fundamental theory, classical thermodynamics appears 
today as a meta-theory applicable to a subset of all dynamical systems subject 
to particular constraints. Since thermodynamics can be abstracted from the 
underlying physics, thermodynamics emerges from the average dynamics of 
certain chaotic systems with many degrees of freedom. 

As consequence, when applicable thermodynamical quantities and their prop- 
erties, in particular entropy, should be derivable from the underlying dynamics. 
Contrary to quantities like time, space, and mass, thermodynamical quantities 
do not need to be assumed to be fundamental, since they can be represented 
on computer by properly averaging the system dynamical quantities. 

1.3 Dynamical systems 

In the following we restrict our considerations on classical systems for which 
the relationship between chaos and integrability on the statistical behaviour is 
much better understood than in quantum systems. Also classical systems can 
be simulated on computers much more easily than quantum systems, which 
helps forging a mental model of their statistical behaviour. Much progress has 
been achieved in the last decades about understanding the correspondence of 
classical chaos, thermodynamics and quantum mechanics (e.g. Gemmer et al. 
2004), but to not complicate the discussion we will restrict it below on classical 
systems. 

A general dynamical system is often defined by the phase space {z}, a set of 
real vectors in a A^-dimensional Euclidean space, generated by a set of first 
order ordinary differential equations 

k=F(z,t). (1) 

The set of all dynamical systems is a too broad set to build thermodynam- 
ics upon. Dissipative systems have attractors, therefore their asymptotic be- 
haviour is restricted to these attractors, which may be as simple as a single 
point in phases space. Richer attractors are limit cycles, and still richer attrac- 
tors, the strange attractors, may have fractal dimension. Reversing time in a 
dissipative system produces an anti-dissipative system which tends asymptot- 
ically to explore the boundaries of the allowed phase space. 

Naturally some systems may have some dissipative phase space coordinates 
and some other anti-dissipative coordinates. Among such systems the Hamil- 
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V = Non Hamiitonian 




Fig. 1. The set of all dynamical systems includes the dissipative (V < 0) and 
anti-dissipative (V > 0) system subsets. Their complement is the subset of volume 
conserving (V = 0) systems, in which the Hamiltonian systems are a subset. Inte- 
grable and ergodic systems are distinct subsets overlapping the previous subsets. 
Classical statistical physics is restricted to the ergodic subset, while deterministic 
physics is restricted to the integrable subset. 

tonian systems form a special subset where the phase space variables may 
be chosen in conjugate pairs in such a way that each of them conserves its 
phase space area. The symplectic structure of phase space follows from the 
Hamiltonian H (p, q) and the particular form of Hamilton's equation of motion. 
Hamiltonian systems still form a too broad set for classical thermodynamics, 
since some of these systems have long range interactions for which energy is 
not scaling proportionally to volume or mass. Others Hamiltonian systems are 
fully integrable, that is decomposable into independent one-degree of freedom 
Hamiltonian systems (consisting each of one pair of conjugate variables). Due 
to the area conservation constraint, such systems cannot be chaotic. 

However, for many Hamiltonian systems classical statistical mechanics can al- 
ready be used without the assumption of extensivity, as discussed by Gross in 
this volume and elsewhere (e.g. Gross 2001), or Padmanabhan for gravitating 
systems (1990). This is the microcanonical description, which however still 
requires the ergodic hypothesis, which in turn requires a finite phase space 
volume. Further, physical systems are modelled over a finite time over which 
a statistical description makes sense. This often means that the strength of 
chaos, measured by the rate of exponential divergence of neighbouring trajec- 
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tories, must be large enough to effectively provides an sufficiently good ergodic 
behaviour over the physically relevant finite time. 

The Fermi-Pasta-Ulam (1955) oscillator model first showed that not any me- 
chanical system is ergodic, and subsequent works on dynamical systems such 
as the Kolmogorov-Arnold-Moser (KAM) theorem provided many more indi- 
cations that most dynamical systems are partly chaotic, or semi-ergodic. Thus 
typical Hamiltonian systems are in general not ergodic, therefore the ergodic 
hypothesis is not necessarily valid in integrable or semi-ergodic systems. 

In integrable systems almost all the orbits are quasi-periodic, therefore the 
initial correlations given by the initial conditions are never erased. Although 
exactly integrable systems are exceptional, fully ergodic systems are also rare. 
Generic Hamiltonian systems possess phase space regions with local quasi- 
integrals, which are nevertheless riddled by a web of resonances and chaos. 
In such semi-ergodic systems the accessible phase space over finite but long 
"quasi-asymptotic" times can be confined to a subset of the non-integrable 
phase space that depends finely on the particular initial conditions. This leads 
to intermittent behaviours. In such systems the usual statistical mechanical 
tools are not appropriate to describe steady states, because sudden large fluc- 
tuations are possible after arbitrary long times. 

In a mechanical description of a physical system the basic properties that 
imprint global symmetries are the invariant functions of the coordinates, the 
global integrals of motion. Since in all time independent mechanical systems 
the total energy function E = H(p(t), q(t)) is the associated invariant, it is not 
astonishing that energy plays a key role in thermodynamics. Any invariant of 
the system introduces a geometrical constraint in phase space that is obeyed by 
the system, despite the presence of many more remaining fully unconstrained 
degrees of freedom. 

Here enters the crucial property that most systems where thermodynamics 
works sufficiently well must have: most of the degrees of freedom must be 
strongly chaotic, which means that beyond a relaxation time the unavoidable 
and not accountable perturbations of the rest of the Universe on the modelled 
system can no longer be neglected. Beyond this time a deterministic descrip- 
tion of the system is impossible. If the perturbations are uncorrelated with the 
system degrees of freedom, then rapidly a strongly chaotic system will have 
an equal probability to visit almost all points of the non-integrable subspace 
of phase space. Therefore the ergodic principle in statistical mechanics only 
applies over the variables belonging to the sufficiently chaotic degrees of free- 
dom. Then over these variables averaging works well. In contrast, each strict 
integral of motion lowers the dimension of the available phase space by one. 

As an example of system with non-thermodynamical behaviour, for most ini- 
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tial conditions a gravitational 10-body system dissolves in a few crossing times, 
while for particular initial conditions like the present solar system ones the 
same 10-body system may stay bound over billions of crossing times. The so- 
lar system is therefore a prime case of semi-ergodic system on which classical 
thermodynamics does not apply. Over time-scale of a few milli on years only a 
small part of the degrees of freedom behave in a chaotic way ([Laskarl . [1994), 
indeed the other degrees of freedom are constraint not only by the 10 global 
classical integrals of motions, but also by additional local constraints that act 
almost like additional integrals of motion over very long finite time-scales. 



A natural global quantity of a bounded mechanical Hamiltonian system fol- 
lowing the orbit j.P(t), Q(t)\ is the total visited phase space volume 



r = hm Y t J dt ' I dN v dN P 

-t 



8(q-Q(t'j)8(p-P(t'] 



(2) 



If the orbit were ergodic, almost all the points at constant value of E of H (p, q) 
would be included in V. This constant E volume can also be expressed as 

T E = Jd N qd N p5(E-H(p,q)) . (3) 



Every additional independent invariant function Ik{p,q) restricts the avail- 
able space space by one dimension in general systems, by two dimensions in 
Hamiltonian systems due to the symplectic constraint. If the number of chaotic 
degrees of freedom exceeds by far the number of invariant functions, the fa- 
miliar situation for gases, the total accessible phase space volume remains a 
good approximation of the total phase space volume. This quantity, taken in 
the log, is, up to a constant, nothing else as Boltzmann's original entropy 
S = fcslnlV. For a statistical mechanics based on a well defined dynamics, 
such an entropy has a well defined sense. Without such a definition based 
on the particular system microdynamics, entropy remains a non-measurable 
quantity with obscure physical meaning. 1 

Explicitely, if the Hamiltonian system has N degrees of freedom, M < N 
global independent integrals Ik(Pi,qi) with values J*,, and if the remaining 
accessible phase space is ergodic, then the log of the accessible volume, 

„ M 

r ergod ic = UHJk- hip, q)) dq^dpf , (4) 
J fc=i 



1 John von Neumann once advised Claude Shannon to call "entropy" the measure 
of information: "No o ne knows what entropy really is, so in a debate you will always 



have the advantage" ( Tribus Irvinel . 197lh . 
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can be used as the system microcanonical entropy, on which a thermodynamics 
can be built. 

1.4 Generalisations 

From there it is straightforward to generalise these considerations for a wider 
class of dynamical systems. First, Hamiltonian systems form a subset of vol- 
ume preserving systems. By taking a Cartesian product of simpler volume 
preserving systems and combining the variables by regular functions one can 
build more general coupled dynamical systems with an odd or even number of 
variables that are still volume preserving, and possess a number of global inte- 
grals of motion. Like for Hamiltonian systems the remaining accessible phase 
space however may or may not be effectively ergodic due to chaos. In the first 
case the entropy can be defined and be used to built upon a microcanonical 
thermodynamics. 

— # — * 

For dissipative dynamical systems z — F(z) with orbit Z(t), the asymptotic 
attractors may still represent a non-trivial invariant subset of phase space 
which is eventually ergodic. For such systems the asymptotic thermodynamics 
is still a well defined procedure, just taken in the asymptotic limit of large time, 



Actually, most steady natural systems can be supposed to follow this route. 
Starting from non-equilibrium states they converge toward asymptotic average 
steady states due to dissipation, and are maintained on their respective final 
attractors by global integrals. 

1.5 The reason for extensive variables 

An important point to explicit is the reason why thermodynamical variables 
in the canonical description are supposed to be either intensive, i.e., position 
invariant, or extensive, i.e., scaling linearly with the spatial volume V. For 
some reason, rarely discussed in depth, the thermodynamical variable Xj of a 
system with many degrees of freedom scales as 



where a = or a = 1. To assume such scaling relations is a crucial step leading 
to the usual characteristics of classical thermodynamics. This symmetry allows 




*l+*2 




ti 



Aj oc V a , 



(6) 
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to describe extensive systems with uniform prescriptions encouraging the belief 
that extensive thermodynamics is universally applicable. 

Let us apply the considerations of previous section to explain why systems with 
a large number of degrees of freedom tend toward uniform density above some 
microscopic scale. The smooth average distribution of particles is necessary to 
describe macroscopic variables as differentiable, which is then required to allow 
the use of differential equations for describing globally out of equilibrium fluids 
only locally in thermal equilibrium. In such systems the particle Lagrangian 
£ for the space coordinates {x{} takes the usual form where T is the kinetic 
energy and U the potential energy, 

C = T&)-U(2i), i = l,...N. (7) 



In systems where the particle interaction energy is negligible with respect to 
their kinetic energy, |T| ^> \U\. The Lagrangian becomes then independent 
on the space coordinate afj. Therefore any integral of the system and subs- 
system should adopt the same space and orientation invariance: the system of 
particles must tend toward a uniform spatial distribution of the particles. In 
velocity space the particles are constraint by a finite total energy, a spherical 
symmetry of the distribution function in velocity space with a finite second 
moment follows; for distinguishable particles the Maxwell-Boltzmann distri- 
bution follows. Therefore the uniform smooth mass density in classical gases 
and fluids in statistical equilibrium follows from the corresponding invariance 
of the Lagrangian and from the chaotic nature of motion of particles. 2 The 
extensivity property follows then from the microscopic dynamical properties 
of chaotic particles subject to a mutual interaction energy small in regards of 
their kinetic energy. The intensive variables, for example the mass density p, 
derive from the ratio of two extensive variables, such as the ratio of the mass 
and the volume. 

At this point it is clear that when the interaction potential energy is not neg- 
ligible, at sufficiently low temperature, the rules of cla ssical thermodynam - 
ics can break down. Indeed, as stressed for example by iLvnden-Bell 



at sufficiently low temperature a phase transition occurs, the micro canonical 
specific heat becomes negative. The system is then thermally unstable, the 
homogeneity invariance breaks down due to coexisting phases. 

At still lower temperature the interaction energy in the Lagrangian dominates; 
a crystal may form, in which case most of the degrees of freedom become 
frozen. The continuous translational invariance is then broken and replaced 



2 Although the perfect gas model neglects formally the interaction term in the 
Lagrangian, collisions are essential in real gas for making the molecular motion 
chaotic. Otherwise the chaos would depend entirely on the boundary conditions. 
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by a discrete lattice invariance. The remaining small motions of an ensemble of 
molecules around the lattice energy minima form again an extensive dynamical 
system which can be described with usual statistical physics at the condition 
that the small motions around the lattice minima are nevertheless chaotic, 
unlike the Fermi-Pasta-Ulam oscillators. 



2 Non-extensive systems 

However many cold systems with strong interaction potential do not develop 
regular lattices, among them glasses, or the gravitational iV-body systems. 
Since the statistical behaviour of gravitational systems is observed to display 
similar properties across galaxies and the Universe, a statistical description of 
them can nevertheless be expected. 

An important additional symmetry in systems of particles such as gravita- 
tional systems is that the long range interaction term is scale invariant. In 
other terms the potential energy follows the symmetry of p-homogeneous func- 
tions, 

U(Xxi) = \ p U{xi) , (8) 

where A is any real positive number, and p is the power law exponent of 
the 2-particle potential; p = — 1 for gravitational systems. For p sufficiently 
negative the potential is short ranged, so we should expect to recover perfect 
gas behaviour at sufficiently negative p. 

So if we follow the above reasoning, a system of strongly interacting particles 
might display statistical behaviour following different symmetries than the 
ones adopted in extensive systems: instead of spatial invariance, systems may 
be characterised by scale-invariant properties. A priori nothing obliges us to 
assume a smooth, differentiable average particle distributions. Distributions 
in the mathematical sense, such as fractal or statistically scale-invariant dis- 
tributions are actually suggested by many physical systems, as popularised by 
Mandelbrot (1982). Such distributions are in general non-differentiable, but 
can still be represented by an ensemble of points, i.e., a superposition of Dirac 
distributions. 

So the problem that we want to address is to describe statistically such systems 
without adopting the usual assumptions of thermodynamics, but using as 
much as possible the known constraints from the microdynamics. Clearly we 
do not want to start with the assumption of extensivity, but rather take the 
symmetries of the Lagrangian as a guide to describe it in a statistical way. 
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Instead of requiring microscopic extensivity which justifies smooth microscopic 
distributions just above the mean-free path scales, we want to allow instead 
scale-invariant hierarchical systems, since such systems appear frequently in 
nature. The consequence is that the usual mathematical tools mostly based 
on differential equations are no longer suited. 

The main general statistical tool of particle dynamics is the virial theorem, 
which applies to any, smooth or non-smooth, distribution of particles. In the 
literature the virial theorem is often derived by averaging the fluid equations of 
motions. This approach assumes from the start that the distribution function 
is regular. However the more general route not requiring regularity is to use 
the Lagrange- Jacobi identity valid for any set of point mass particles. 

The polar semi-moment of inertia / about the origin of a system of particles 
of mass rrii 

J=i]Tm^, (9) 

i 



characterises the system average mass extension in space. Differentiating / 
with respect to time, the acceleration of this global scalar quantity is seen to 
be linked to the system kinetic energy and the total moment of radial forces, 
where F { = rriiXi, 

i = ^2m i x i +^2F i -x i . (10) 

i i 



If the forces derive from a sum of p^-homogeneous potentials U k , 

'x = -J2VU k (x l ), where U k (Xx) = \ Pk U k (x) , (11) 
k 



then the total moment of force is proportional to a sum of potential energies 
weighted with the homogeneity indexes p k , thus 

I = 2T(x 2 i )-J2PkU k (x' i ). (12) 



This is the Lagrange- Jacobi identity, slightly generalised with a sum of homo- 
geneous potentials. 

If the system is confined by a pressure at the boundary, this pressure has 
the same effect that a uniform negative kinetic energy distributed inside the 
system volume V. Since a non-relativistic kinetic pressure amounts to 2/3 of 
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13 

Fig. 2. A hierarchical 2 4 -body system, 
the kinetic energy density, we can write in virial equilibrium (/ = 0) that 

2T{x1) -J2PkU k {xi) = ZPcxtV . (13) 
k 

3 Hierarchical systems 

3.1 The hierarchical N L -body problem 

The gravitational TV-body problem has played a central role in the develop- 
ment of classical mechanics and the theory of dynamical systems. It was seen 
for over 2 centuries as a prime example of the deterministic physics devel- 
oped by Newton and followers. The main reason is, first, that the 2-body case 
is fully integrable, so computable over an arbitrary long time, and, second, 
that the planets move on orbits that are perturbations of the 2-body problem. 
Therefore the solar system remains well computable over centuries. But the 
works on the 3-body problem by Brun and Poincare in the late 19 th century 
showed that the general iV-body problem was of a different kind, not so eas- 
ily computable. The notion of exponential divergence of nearby trajectories 
leading to the modern notion of chaos came then into existence. 

Hierarchical iV-body problems have been rarely studied (e.g., Nugeyre & Bou- 
vier 1981). Yet such systems may provide others configurations than the plan- 
etary systems that are close to be integrable. Let us consider a bound 2-body 
gravitational system. If now we replace each body by an identical but scaled 
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down 2-body problem at its centre of mass, and then recursively over a number 
L of levels, then we have built a hierarchical 2 L -body problem (Fig. 2). This 
system is close to be integrable, since each 2-body subsystem is perturbed 
by the inner and the outer other 2-body systems. In order to avoid possible 
resonances we have the freedom to choose the scale ratio is such a way as to 
avoid rational ratios of orbital periods. Then in the limit of a large scale ratio 
the respective perturbations by the other bodies on each two body system is 
small, eventually negligible for at least the finite time we wish to consider such 
a system. 

3.2 Geometrical relations 

Let us first enumerate some simple but little known geometrical properties of 
mass particles arranged in a hierarchical way (Pfenniger & Combes 1994). 

As first hypothesis we consider a system of mass mi made of N clumps of 
mass mi_i, made themselves of N sub-clumps m L _2, and so on over a number 
of levels. At any level L > 0, 



In a regular hierarchical system the ratio of masses between adjacent levels is 
supposed to scale with spatial size according to a power law d, 



mi-i \ r L-iJ 

Since mass is positive, its cumulative value within a sphere centred on a mass 
point must monotonously increase with the sphere radius, so d > 0. In a 
3-dimensional space, the positive cumulative mass cannot grow faster than 
volume, therefore d < 3. 

As consequence there is a direct relationship between the adjacent level size 
ratio, N and d, 



m L = Ntjil 



-i 



(14) 




(15) 



= N 1/d . 



(16) 



TL-l 



It follows straightforwardly that the average mass density scales as 




(17) 
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which means that the average density increases at small scale for all permitted 
d < 3. 

In contrast, the projected mass density T, L scales as 
S L / m L /r 2 \ ( r L \ d ~ 2 

e^T = l^wfcj = {—J ' (18) 

Thus the critical exponent d = 2 separates the projection of fractal objects 
of different d (see also Falconer 1990). The immediate consequence for astro- 
nomical observations restricted by projected views of cosmic structures is that 
if a fractal object has 3 > d > 2 then the sky is well filled by it, like common 
fog does. However, when d < 2 most of the mass is projected in small spots, 
leaving most of the sky empty. Therefore a fractal distribution of galaxies with 
d < 2 (Coleman & Pietronero 1992) is a possible (partial) explanation of the 
blackness of the Universe, the famous de Cheseaux-Olbers paradox. 



3.3 Dynamical relations 



Hierarchical distributions of gas are observed in the interstellar medium with- 
out a clear understanding of the overall physics. It seems clear that in a galaxy 
gravity plays a non-negligible role at most scales. Typically the interstellar gas 
fills the volume of a galaxy in a very inhomogeneous way well described by 
fractal geometry. The matter is fragmented in clumps themselves composed 
of sub-clumps over several order of magnitudes in scale. At each scale of a 
galaxy mass condensations are rather close to virial equilibrium from sizes 
comparable to the whole galaxy down to stars, passing by various interstellar 
clouds of intermediate sizes. Between the largest galaxy size and the smallest 
stellar size, a range of order 10 12 is covered. 

Generalising the gravitating case to other power law potentials, we adopt the 
potential form between mass particles, 

$( r ) = ^ r p ) (19) 



with p 7^ 0. The limit p — > is well behaved and provides a logarithmic 
potential, that we do not discuss separately. For attracting interactions the 
coupling constant G is positive. The gravitational case occurs for p — — 1. 

Suggested by the system hierarchical organisation, we approximate the poten- 
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Fig. 3. The mass hierarchical model over 4 levels with N = 3. 



tial energy Ul at level L by retaining the first term of a multipolar expansion, 



The potential energy at a given level consists of the summed sub-clumps po- 
tential energies and the potential energy due to the interaction between the 
sub-clumps. The term containing a represents the first correction taking into 
account tidal interactions. Since the scale ratio is constant, the term in paren- 
theses is constant and can be absorbed in a new coupling constant G' . We will 
see that the main results do not depend on the value of G (indeed scale invari- 
ance effects depend on the power law exponent p). For parameters suited to 
the interstellar medium [D m 1 — 2, iV = 5 — 8, Scalo 1990) the approximation 
(20) has been checked to be accurate at the percent level. 

The third hypothesis is of statistical character. Thermodynamics can not be 
used since it excludes long-range interactions. The only remaining statisti- 
cal tool is the virial theorem, more precisely the Lagrange- Jacobi indentity. 
If a scale invariant system finds a statistical equilibrium at each level, then 
the semi-moment of inertia acceleration I of different sub-clumps vanishes in 
average, such that at each level, 



where T L is the total kinetic energy summed up level L, Pl+i is the outer 
kinetic pressure, and Vi ~ r\ the clump volume. Here the outer pressure is 




(20) 



2T L -pU L = 3P L+1 V L , 



(21) 
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purely kinetic: it is given by 2/3 of the kinetic energy density outside the 
sub-clumps. Approximately, 

Pl ~ 3V L -NV L ^ (22) 



which takes into consideration that most of the time distinct clumps do not 
overlap. 



3.4 Explicit solution 



The above recurrences for M L , V L , U L and T L in Eqs. (14,15,20,21,22) can be 
solved with little algebra exactly in finite terms, as detailed below. 

First we define a few convenient terms, and state simple constraints. The 
lowest level virial ratio, 

(23, 



is obviously constraint from Eq. (21) divided by 2T L in the interval < (3 < 1, 
because Tl, pUi, Vl, and Pl are all positive. When j3 — > the lowest level is 
confined by outer pressure, when f3 — > 1, by internal attraction. 

The lowest level energy is E = T + U . The marginally bound state occurs 
when E = 0. Such a state is possible only when 

< -p/2 < 1 . (24) 



The scale ratios 

x = — , and X = x L , (25) 
r L -i 



are by definition > 1, since as shown in Eq. (16), x = N l l d > 1. 
In terms of x and X we have, 

^ = X d , and ^ = X 3 . (26) 
m V 

It just means that the mass is proportional to the power d of the size, while 
the volume to the power 3. Incidentally, the factor X d ~ 3 is the density ratio 
Pl /po- 
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Eq. (20) is a linear recurrence relation for Ul that has a simple solution: 

U L /m L (Xx) d+P - 1 



U /m ~ x d+ P - 1 ' 



L > 0. 



(27) 



Some more effort is required to solve Eq. (21) by substitutions and induction. 
We obtain, 



T L /m L 
T /m 



x d+p_l 



{Xxf^-{Xx) d+p x 3 ^-X d+p 



x dJrP 



1 1 

— d ,gd-\-p 



(28) 



The equivalent of the kinetic pressure at each level is proportional to the 
square of the velocity dispersion vl given by, 



-m L v L = T L , 



= T T , — NT, 



L-l 



L = 0, 
L > 1 , 



(29) 
(30) 



which means that the kinetic pressure is the difference between two adjacent 
levels of the summed kinetic energies down to the lowest level. Replacing with 
Eq.(28), we obtain, 



(x 3 -" - l) 



-) 

x J 



3-d 



P 



X 



3-d 



X d+ P 



rg3 — d rgd+p 



L > 1 



(31) 



3.4-1 Physical domain of p and d 

The above hierarchical model makes sense in the intervals < d < 3, for 
geometrical reason, and p < if we wish to consider interactions of attractive 
type not growing too fast with distance. As already understood by Newton, 
even the gravitating case p — — 1 leads to conceptual problems when d = 3 
(the force at any point in a uniform infinite medium is a sum of cancelling 
diverging forces). 

However, even in this restricted domain not any combination of p and d is 
physically feasible. One has to ensure that the positive quantities are indeed 
positive, and do not diverge to excessive values. 

First one can easily check in Eq.[27] that U L /U does not change sign because 
X > x > 1, and at d + p = the behaviour is regular. 

The kinetic energy and the squared velocity (or the pressure) should remain 
positive and finite. In fact, by looking at Eq.[30] ones sees that it suffices to 
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N,|M 





d 




i.) 
-1 
p-2 
-3 
-4 



N= 10000 



d 



Fig. 4. The ratio v\/vq, for I = 12, in the plane (d,p) takes either very large values 
exceeding 10 4 in the shaded regions at low d, or negative values in the white region 
at large d, or physical values in the narrow intermediate region, for different N. In 
all cases the low level boundary condition is (3 = min(l, —p/2), which produces the 
discontinuity at p = — 2. 



require the positivity of the squared velocity at each level to ensure that the 
summed kinetic energy is also positive. 

A denominator vanishes in Eq. [28] when d+p —>■ 0. But more important, when 
3 — d — > d + p other denominators in Eqs.[28],[31] vanish with corresponding 
numerators containing cancelling powers of X. This latter singularity is cru- 
cial because it means that as the number of levels of the system increases, a 
shrinking domain of parameter space is able to correspond to accessible physi- 
cal states. As shown in Fig. (4), on one side of the singularity v\ is negative and 
on the other side it takes extremely large values, such that the only physical 
sensible region is within the singularity as L — > oo. 
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Since v\ can change its sign, we obtain a limit of physically admissible cases. 
Solving Eq.[31] for f3 we obtain a constraint for f3, noting a = (d+p) — (3 — d), 



v 2 L >0 =>- p< 



l-x a 
1-X a 



(32) 



Since X > x > 1, the critical value /3 cr i t is in the interval (0 — 1) only when 
a > 0. To not be subject to this restriction on (3, d must obey 



For hierarchical gravitating systems with d > 2, a range of (3 < 1 is not 
allowed. At the lowest level an outer pressure is required. 

3.5 Large L limit 

In ordinary systems the thermodynamical limit means that the boundary con- 
ditions, the "box", takes an infinite dimension, V — > oo and M — > oo, but 
M/V — > const < oo, or the number of particles N tends toward infinity at con- 
stant average density. The mass and volume per particle are kept constant. In 
other cases, particularly relevant in gravitating systems such as stars in which 
local thermodynamics still works, the total mass and volume is taken as fixed, 
but the number of particles is increased to infinity by decreasing the mass m 
and volume v per particle to zero, keeping m/v constant. 

In hierarchical systems, one can also increase the number of particles, or de- 
grees of freedom, toward infinity at either the large or small scales, by in- 
creasing the number of levels L. The boundary conditions, fixing typically the 
temperature in a box, can be replaced by the conditions at either small or 
large scale depending on the problem. 

Here we take the analogue of a "thermodynamical limit" by L — > oo, or 
X = x L — > oo. Then the shrinking of the range of physical solutions compatible 
with the virial theorem leads to the constraint, 




(33) 



d = 



3 — p 



(34) 



ln(l-/3) " 
In N 



Equivalently, one can use 



(3 = 1 - x a , 



a = 



(d + p) - {3-d) . 



(35) 
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For example for p = — 1, d < 2 in any case. For systems not too confined 
by the outer pressure (1/2 < (3 < 1) and iV > 5, we have d < 1.7, while 
systems with more outer pressure ("pressure confined clouds") increase d up 
to 2. Therefore, we conclude that hierarchical gravitating systems in statistical 
equilibrium can indeed exist, but with d < 2. 

3.6 Velocity scaling 

In the large L limit, when using Eq. [35] the scale- velocity dispersion in Eq.[31] 
takes the exact scale-free form 

f=(l-^)^. ^1- <*» 

In other words 

K i d + p 
v oc r , where k = . (37) 

For fractal ISM cold clouds with substantial ambient pressure, (3 1/2, and 
with iV = 5 — 10 (Scalo 1990), we expect d 1.7 and « 0.35, which is 
compatible with Larson's (1981) size-linewidth relationship. 



4 Conclusions 

Once one realises that despite its enormous success, classical statistical me- 
chanics is not universal, but mostly applicable to extensive strongly chaotic 
systems, a huge field of investigations opens up. Classes of systems, including 
formal or non-physical systems, may eventually be found suited for a statis- 
tical description of their steady states provided their dynamics is sufficiently 
chaotic to allow the use of an ergodic hypothesis. However, many semi-ergodic 
systems exist with intermittent behaviour for which there is little motivation 
to expect that a statistical description of steady states is appropriate. 

The non-extensive gravitational systems are the rule in astrophysics, and many 
observational evidences exist that suggest that gravitational systems may fol- 
low statistical rules. For example laws such as the r 1 / 4 radial density profile of 
elliptical galaxies, or the exponential profiles of stellar galactic disks are well 
known, and reproducible in computer simulations, but not fully explained by 
a statistical model. 
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Here as a illustrating example we have developed a statistical description 
of N L hierarchical systems subject to power law interaction without using 
concepts of classical statistical mechanics. Using the virial theorem instead 
one can already find new constraints that emerge for a large number of bodies. 
As application, Larson's empirical velocity scaling relationship in molecular 
clouds is retrieved. 
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